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We study elastic systems such as interfaces or lattices pinned by correlated quenched disorder 
considering two different types of correlations: generalized columnar disorder and quenched defects 
correlated as ~ a;"" for large separation x. Using functional renormalization group methods, we 
obtain the critical exponents to two-loop order and calculate the response to a transverse field h. 
The correlated disorder violates the statistical tilt symmetry resulting in nonlinear response to a 
tilt. Elastic systems with columnar disorder exhibit a transverse Meissner effect: disorder generates 
the critical field he below which there is no response to a tilt and above which the tilt angle behaves 
as 1? ~ (ft — hc)"^ with a universal exponent (/> < 1. This describes the destruction of a weak Bose 
glass in type-Il superconductors with columnar disorder caused by tilt of the magnetic field. For 
isotropic long-range correlated disorder, the linear tilt modulus vanishes at small fields leading to a 
power- law response i9 ~ /i'* with > 1. The obtained results are applied to the Kardar-Parisi-Zhang 
equation with temporally correlated noise. 

PACS numbers: 74.25.Qt, 75.60.Ch, 81.10.-h 



I. INTRODUCTION 

Elastic objects in disordered media are a fruitful con- 
cept to study diverse physical systems such as domain 
walls in ferromagnets,— charge density waves in solids 
(CDW)^ and vortices in type-II superconductors^ In 
all these systems, the interplay between elasticity, which 
tends to keep the object ordered (flat or periodic), and 
disorder, which induces distortions, produces a compli- 
cated energy landscape i^i^ This leads to rich glassy be- 
havior. For instance, at low temperature, weak defects in 
a crystal of type-II superconductor, such as oxygen va- 
cancies, can collectively pin the flux lines in the so-called 
Bragg glass statei^ Vortex pinning prevents the dissipa- 
tion of energy, and thus, its understanding has a great 
importance for applications. It was observed in exper- 
iments that columnar defects produced in the underly- 
ing lattice of superconductors by heavy ion irradiation 
can significantly enhance vortex pinning.^ Nelson and 
Vinokur^ mapped the problem of flux lines pinned by 
columnar defects onto the quantum problem of bosons 
with uncorrelated quenched disorder in one dimension 
less. The mapping predicts a low temperature "strong" 
Bose-glass phase which corresponds to the localization of 
bosons in a random potential provided the longitudinal 
applied field TJy is weak enough to create vortices with 
density smaller than the density of pins. For larger H\\ , 
the Bose-glass can coexist with a resistive liquid of inter- 
stitial vortices which, it is argued, can freeze upon cool- 
ing into a collectively pinned weak Bose-glass phase. At 
low tilts of the applied magnetic field relative to the par- 
allel columnar defects, flux lines remain localized along 
the defects, so that vortices are characterized by an in- 
finite tilt modulus. This phenomenon which is known 
as the transverse Meissner effect has been extensively 
studied experimentallyf^'^l Vortices undergo a delocaliza- 
tion transition to a flux liquid state at some finite criti- 
cal mismatch angle i^c between the applied field and the 



direction of defect alignment, i.e., at some finite trans- 
verse field HI. The schematic phase diagram is shown 
in Fig. [1] The breakdown of the transverse Meissner ef- 
fect above can be described by 

^ {H^ - Hl)t (1) 

where B± is the transverse magnetic induction due to the 
tilted flux lines. Heuristic arguments of Ref. JJ, based on 
kink statistics predict (jj = 1/2 in d = 1 + I dimensions 
and (/) = 3/2inrf = 2+ l. However, experiments on 
a bulk superconductor (d — 3) with columnar disorder 
find <j> ~ 0.5,^^ while the strong-randomness real-space 
renormalization group suggests cj) = 1 in d = that 
is in disagreement with the predictions based on kink 
statistics. Thus, further investigations are needed. 

The theoretical advances for elastic objects in dis- 
ordered media are achieved by developing two general 
methods: the Gaussian variational approximation (GVA) 
and the functional renormalization group (FRG). GVA 
relies on the replica method allowing for the replica sym- 
metry breakingii^ It is exact in the mean field limit, i.e., 
in the limit of a large number of components. FRG is 
a perturbative renormalization group method which is 
able to handle infinite number of relevant operator.— 
Simple scaling arguments show that the large-scale prop- 
erties of a c?-dimensional elastic system are governed by 
uncorrelated disorder in d < due = 4. In particular, 
displacements grow unboundedly with distance, result- 
ing in a roughness of interfaces or distortions of peri- 
odic structures. The problem is notably difficult due to 
the so-called dimensional reduction which states that a 
d-dimensional disordered system at zero temperature is 
equivalent to all orders in perturbation theory to a pure 
system in d — 2 dimensions at finite temperature. How- 
ever, metastability renders the zero-temperature pertur- 
bation theory useless: it breaks down on scales larger 
than the so-called Larkin lengthfil The peculiarity of 
the problem is that for d < due there is an infinite 
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FIG. 1: Phase diagram of flux lines with columnar disorder at 
a fixed applied field _ff|| parallel to the columns as a function 
of transverse field H± and temperature T. The Bose glass 
corresponds to vortex localization and exhibits a transverse 
Meissner effect. The transition to liquid phase occurs at the 
critical field H1{T). In (1+1) dimensions, the temperatures 
Tbg and T* are related to the special values of the Luttinger 
interaction parameter: g{TBG) ~ 3/2 and g{T*) = 1. The 
inset shows the response to a transverse field in the vicinity 
of transition. 



set of relevant operators. They can be parametrized 
by a function which is nothing but the disorder cor- 
relator. The renormalized disorder correlator becomes 
a nonanalytic function beyond the Larkin scale. The 
appearance of a nonanalyticity in the form of a cusp 
at the origin is related to metastability, and nicely ac- 
counts for the generation of a threshold force at the de- 
pinning transition . ^^'^^1^°'^-'^ It was recently shown that 
FRG can unambiguously be extended to higher loop or- 
der so that the underlining nonanalytic field theory is 
probably renormalizable to all orders . ^^'^'^ '^ Although 
the two methods, GVA and FRG, are very different, they 
provide a fairly consistent picture and recently a rela- 
tion between them was established.^^ There is also good 
agreement with results of numerical simulations, not only 
for critical exponent o^^'^^i^^ but also for distributions of 
observables22i^ and the effective actionj2i 

The FRG techniques were also applied to pinning of 
elastic systems by columnar disorder j^i^^i^ The mod- 
els studied by FRG, though that may be more directly 
applicable to systems such as charge density waves or do- 
main walls, exhibit many features of the Bose-glass phase 
of type-II superconductors. In particular, they demon- 
strate the absence of a response to a weak transverse field 
and provide a way to compute the exponent cf).^'^ How- 
ever, since FRG intrinsically assumes collective pinning, 
it also predicts a slow algebraic decay of translational or- 
der, that is not expected in the strong Bose-glass state 
when each vortex is pinned by a single columnar pin. 
Thus, the FRG is able to handle only a weak Bose-glass 
phase, exhibiting both the transverse Meissner effects 
and the Bragg peaks. 

In the present paper, we extend the FRG studies to 
two-loop order. We also extend to two-loop order our re- 
cent work^ on the elastic objects in the presence of long- 
range (LR) correlated disorder with correlations decaying 
with distance as a power law. This type of disorder can 



be induced, for example, by the presence of extended de- 
fects with random orientations. In particular, we address 
the question of the response to a tilting field and compare 
the effects produced by different types of disorder corre- 
lations. The outline of this paper is as follows. Section HIl 
introduces the models of elastic objects in the presence 
of generalized columnar and LR-correlated disorder. In 
Sec. mil we study the model with LR correlated disorder 
using FRG up to two-loop order. In Sec. IIVI we con- 
sider the response of elastic objects to a tilting field and 
discuss the relation to the quantum problem of interact- 
ing disordered bosons. In Sec. |Vl we revise the problem 
of surface growth with temporally correlated noise using 
the results obtained in the previous sections. 



II. MODELS WITH CORRELATED DISORDER 

The configuration of elastic object embedded in a 
D-dimensional space can be parametrized by an iV- 
component displacement field Ux, where x belongs to 
the d-dimensional internal space. For instance, a d- 
dimensional domain wall corresponds to d = D — 1 and 

= 1, vortices in a bulk superconducter to d = D = 3 
and N — 2, and vortices confined in a slab to d = D = 2 
and A^ = 1. In this paper, we restrict our study to the 
case A^ = 1 and elastic objects with short-range elastic- 
ity. In the presence of disorder, the equilibrium behavior 
of the elastic object is defined by the Hamiltonian 



n\u] 



-{Vuxf + V{x,Ux) 



(2) 



where c is the elasticity and V{x^u) is a random Gaus- 
sian potential, with zero mean and variance that will be 

defined below. We denote everywhere below = / 

and / = / d'^x. The short-scale UV cutoff is im- 
plied at q ~ A and the system size is L. The ran- 
dom potential causes the interface to wander and become 
rough with displacements growing with the distance x as 
C{x) ^ x^^. Here, C is the roughness exponent. Elastic 
periodic structures lose their strict translational order 
and exhibit a slow logarithmic growth of displacements, 
C{x) = Ad^n\x\. Although most results of the paper 
concern the statics at equilibrium, it is instructive to give 
a dynamic formulation of the problem. The driven dy- 
namics of the elastic object in a disordered medium at 
zero temperature can be described by the following over- 
damped equation of motion 



vdtUxt = cV'^Uxt + F{x, Uxt) + /■ 



(3) 



Here, r] is the friction coefficient, F — —duV{x,u) the 
pinning force, and / the applied force. The system un- 
dergoes the so-called depinning transition at the criti- 
cal force /c, which separates sliding and pinned states. 
Upon approaching the depinning transition from the 
sliding state f ^ fc the center-of-mass velocity v = 



3 



^ /x ^t^Kt vanishes as a power law 



(4) 



In the present work, we consider model ([2]) with two dif- 
ferent types of correlated disorder, which are described 
in two subsequent sections. 



A. Generalized columnar disorder 

Real systems often contain extended defects in the 
form of linear dislocations, planar grain boundaries, 
three-dimensional cavities, etc. We consider the model 
with extended defects which can be viewed as a gen- 
eralization of columnar disorder. The defects are e^- 
dimensional objects (hyperplanes) extending through- 
out the whole system along the coordinate X|| and ran- 
domly distributed in the transverse directions x± with 
the concentration taken to be well below the percolation 
limit i^i^i^ The corresponding correlator of the disorder 
potential can be written as 

V{x,u)V{x',u') = R{u~u')6'^-''''{x^-x'j^). (5) 

The case of uncorrelated pointlike disorder corresponds 
to Erf = and the columnar disorder to = 1. For inter- 
faces, one has to distinguish two universality classes: ran- 
dom bond (RB) disorder described by a short-range func- 
tion R{u) and random field (RF) disorder corresponding 
to a function which behaves as R{u) ~ |u| at large u. 
Random periodic (RP) universality class corresponding 
to a periodic function R(u) describes systems such as 
CDW or vortices in d = 1 -I- 1 dimensions. - 

The standard way to average over disorder is the 
replica trick. Introducing n replicas of the original sys- 
tem we derive the replicated Hamiltonian as follows: 



T 



(6) 



where we have added a small mass m providing an in- 
frared cutoff. Replica indices a and b run from 1 to n 
and the properties of the original disordered system can 
be restored in the limit rt — > 0. We explicitly show in 
Hamiltonian ^ that one has to distinguish the longi- 
tudinal and transverse elasticity modules. Even if the 
bare elasticity tensor is isotropic, the effective elasticity 
may not due to the renormalization by anisotropically 
distributed disorder. 



B. Long-range correlated disorder 

In the case of isotropically distributed disorder, the 
power-law correlation is the simplest assumption with 



possibility for scaling behavior with new fixed points 
(FPs) and new critical exponents. The bulk critical 
behavior of systems with RB and RF disorder which 
correlations decay as a power-law x""' was studied in 
Refs. [39,,40...4L42i. The power-law correlation of disor- 
der in the d-dimensional space with exponent a = d — Sd 
can be ascribed to Sd dimensional extended defects ran- 
domly distributed with random orientation. For in- 
stance, a — d corresponds to uncorrelated pointlike de- 
fects, and a — d — 1 {a = d — 2) describes infinite lines 
(planes) of defects with random orientation. The power- 
law correlation with a noninteger value a = d — dj can 
be found in the systems containing fractal-like structures 
with the fractal dimension Here we consider the 

model with LR-correlated disorder introduced in Ref. [s^ 
which is defined by the following disorder correlator: 



V{x,u)V{x',u') ^ Ri{u-u')S'^{x-x') 

+ R2iu~u')g{x~x'), (7) 

with g{x) x~'^. We fix the constant in the Fourier space 
taking g{q) = q°-~'^. The first term in Eq. ([7]) corresponds 
to pointlike disorder with short-range (SR) correlations 
and the second term to LR-correlated disorder. A priori 
we are interested in the case a < d when the correlations 
decay sufficiently slowly, otherwise the disorder is simply 
SR correlated. 

Using the replica trick, we obtain the replicated Hamil- 
tonian 7in[u] and the corresponding action S[u\: 



S[u] 



T 2T 



„ Jx 



-^E / R^i<~<) 

-^E/ ,^2K-4')ff(^-^')- (8) 

a.b •'^2;' 

One could start with model ([5|, setting Ri{u) = 0. How- 
ever, as was shown in Ref. ISSl . a nonzero Ri{u) is gen- 
erated under coarse graining along the FRG flow. Note 
that the functions Ri{u) can themselves be SR, LR, or 
RP. The generalization of these universality classes to 
LR-correlated disorder is discussed in Ref. [S^- 

In the case of uncorrelated disorder, the system ^ ex- 
hibits the so-called statistical tilt symmetry (STS), i.e., 
invariance under transformation Ux — *■ Ux+fx with an ar- 
bitrary function fx- The STS issues that the one-replica 
part of the replicated action, i.e., the elasticity, does not 
get corrected by disorder to all orders. The presence of 
LR-correlated disorder or extended defects destroys the 
STS, and thus allows for the renormalization of elasticity. 

For a non-Gaussian distribution of disorder, higher or- 
der (p > 2) cumulants would generate additional terms 
in the action with factors of l/T^ and free sums over p 
replicas. These terms are irrelevant in the RG sense that 
can be seen by power counting, and thus will be neglected 
from the beginning. 
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FIG. 2: Two-replica one-loop diagrams correcting disorder. 
The dot line corresponds to either SR disorder vertex (dashed 
line) or to LR disorder vertex (wavy line) . Diagrams of types 
(i) and (ii) contribute to SR disorder. Only diagrams of type 
(iii) give corrections to LR disorder. 



III. RENORMALIZATION OF THE MODEL 
WITH LONG-RANGE CORRELATED 
DISORDER 

A. Perturbation theory and diagrammatics 

We now study the scaling behavior of model ([5]) start- 
ing with simple power counting. The elastic term in ac- 



tion ([8]) is invariant under x — > xb, u 
provided T b^-^T with Qt = d-2-H2C-'(/'- Since Qt is 
positive near d = 4 the temperature T is formally irrele- 
vant. The STS would fix "0 = 0, however, this is not the 
case here. ^ and '0 are for now undetermined and their 
actual values will be fixed by the disorder correlators at 
the stable FP. Under the rescaling transformation the dis- 
order correlation functions R\ and i?2 go up by factors 

tively. Thus, in the vicinity of Gaussian FP (Ri = 0), SR 
disorder becomes relevant for C — -0/2 < (4 — d)/4 and 
LR disorder is naively relevant for (' — ■0/2 < (4 — a)/4. A 
posteriori these inequalities are satisfied at the RB and 
RP FPs. For RF disorder, however, power counting sug- 
gests that SR disorder is relevant for C — "0 < (4 — d)/2, 
while LR disorder is relevant for C — < (4 — a)/2M. 

Let us consider the perturbation theory in disorder and 
its diagrammatic representation. In momentum space, 
the quadratic part of action ([5]) gives rise to the free prop- 
agator (■u°'u^,)o = (27r)'^(5'*(q + q')T5abC{q) represented 
graphically by a line: 



= TC{q)5ab = 



T^ab 



(9) 



We will distinguish two different interactions, SR and LR, 
for which we adopt the following splitted diagrammatic 
representation: 



'T 

I 



E 

ab 

E 

ab 



Rii<-0 

R2{U%-UI,) 

2T2 



g{x - x ). 



(10) 



(11) 



Following the standard field theory renormalization pro- 
gram, we compute the effective action and determine 
counter-terms to render the theory UV finite as d, a — s- 4. 
To regularize integrals, we use a generalized dimensional 



regularization with a double expansion in £ = 4 — d 
and 5 = A — a. The effective action T[u] is defined by 
the Legendre transform T[u] = Ju — W[J], >V[J] = u 
of the generating functional for connected correlators 
yV[J] — lnZ[J]. The replicated partition function Z in 
the presence of sources J is given by 



z[j] = I Yl^n- (-'^M + / E 



(12) 



The effective action T[u\ is by definition a generating 
functional of one-particle irreducible vertex functions. 
However, it turns out to be nonanalytic in some direc- 
tions, and therefore, the relying on the expansion in u 
is danger. To overcome these difficulties, we employ the 
formalism of functional diagrams introduced in Ref. [2^ . 
Since the temperature is formally irrelevant, we compute 
the correction to the effective action at T = 0. Analyzing 
UV divergences of the functional diagrams contributing 
to the effective action, we find that the disorder is cor- 
rected only by local parts of two-replica diagrams and 
the elasticity only by one-replica diagrams. 



B. 



Correction to disorder and P functions 



To one-loop order at T = 0, the correction to disorder 
is given by the local parts of the two-replica diagrams 
shown in Fig. [2l The corresponding expressions read 



RUu)' ^ R'(o{u)R'{o{0) 



h 



2 

+ [R'MR'M - R'MR'^om h 
+ \R2o{^fh, 



(13) 



5^R2{u) = -R%{u)R%{0)h^R%{u)R'i,{0)h, (14) 

where we have included factor of 1/cq in RiQ(u). In this 
section, bare parameters are denoted by the subscript 
"0". The one- loop integrals /i, I2 and I2 diverge loga- 
rithmically and for £, 5 — > are given by 



h = 

h = 
h = 



1 



(g2 



„a — d 



(<Z2 



= K4 



2(a-d) 



{q^ + m?Y 



26^ 



-Oil), 
0(1), 
Oil) 



(15) 
(16) 
(17) 



where we have set m = m/y/cQ and Kd is the area of a 
d-dimensional sphere divided by (27r)''. Let us define the 
renormalized dimensionlcss disorder i?, as 



m^Riiu) = Rioiu) + S^Ri{u), 
m^R2{u) = i?2o(u) + 5^R2{u). 



(18) 
(19) 



Note that to one-loop order, there is no correction due 
to the renormalization of elasticity (see below). The (3 



5 



functions are defined as the derivative of Ri{u) with re- 
spect to the mass m at fixed bare disorder RiQ{u). It is 
convenient to rescale the field u by m'' and write the /3 
functions for the function Ri = K 4m~'^^ Ri{um'^) . Drop- 
ping the tilde subscript, the flow equations to one- loop 
order read 






diRi{u) = {e~4:C)Riiu) + CuR[{u) 



diR2{u) 



R'i{u)f + AR'liu), (20) 
[S - 4C)i?2(w) + (uR'^iu) + AR'^iu), (21) 



where A = -[i?i(0) + i?^'(0)] and di := -m^. The 
FPs of flow equations (pO|) and pTjl characterizing differ- 
ent universality classes have been computed numerically 
in Ref. [s^ and the corresponding critical exponents have 
been derived to first order in e and S. The remarkable 
property of the FRG flow is that the LR part of disor- 
der correlator R2(u) remains an analytic function along 
the flow for all universality classes. We will show below 
that due to this feature, one can obtain the critical ex- 
ponents to two-loop order just computing the two-loop 
correction to elasticity and avoiding exhaustive two-loop 
calculations. 



C. Correction to elasticity 

The STS violation causes a renormalization of elastic- 
ity. The first order correction to the single-replica part 
of effective action is expressed by the following diagram: 



2T -X ), 

a 

(22) 

where the bare correlation function C{x) is given by 
Eq. ([9]). Using the short distance expansion 



d 

E( 

i=l 



4) 



dxi 



(23) 



and identifying the terms of the kind — (Vu°)'^/2T as a 
correction to elasticity, we find 

= ^Mo'(O) / x^9{x)C{x) 



(4).„.<y- 



coi?ro'(o) 



45 



—m 



+ 0(6,6), (24) 



where in the last line we have included K4 /cq in a redefi- 
nition of i?2o(w). Since Eq. is finite for e,S 0, the 
elasticity does not get corrected to one-loop order. 

We now turn to the two-loop corrections. The three 
different sets of diagrams contributing to elasticity are 
depicted in Fig. [31 The details of calculations are given in 
the Appendix. Summing up all contributions, we arrive 



FIG. 3: Single-replica two-loop diagrams correcting elastic- 
ity. The dotted line corresponds to either SR disorder vertex 
(a, /3 = 1) or to LR disorder vertex (a, /3 — 2). The solid line 
corresponds to bare correlation functions C(y). The corre- 
sponding expressions are computed in the Appendix. 



at 



Co 



= i?'/o(0)4?(0) 



S 1 



^^2o(0)4o'(0)™ 



^-2(5 



e + S ie 

e-6 



8^2 



0(1). (25) 



To render the poles in e and 5, we introduce the renor- 
malization group Z factor as follows: 



C^Zc[Rl,R2] ^Cq. 

The exponent ip is given then by 



V' = -m—\nZJRi,R2] 
am 



(26) 



(27) 



where subscript indicates a derivative at constant bare 
parameters. Taking the derivative with respect to the 
mass, we obtain 



-m- — In Zc 
dm 



^^Ri'^m3-e)m~' 



+R'l,{0)R^i^iO)m 
+^2o(0)4oH0)m 



-25 



.5 



4e 

e 



45 



(28) 



To calculate V', we have to express the bare disorder via 
renormalized one as follows: 



Mo^(O) 



M^)(0)-fi?'/(0)i?f (0)1 
i?2(0)i?f (0)i 



(29) 



Substituting Eq. in Eq. ([281), we find that the leading 
two-loop corrections are exactly canceled by the counter- 
terms, so that we leave with 



1 



^ = _l(<5_£)i?f(0). 



(30) 



The finite part of the single-replica two-loop diagrams 
(1^5)1 is expected to correct elasticity at three-loop order. 
Hence, we argue that the perturbation theory for this 
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FIG. 4: (Color online) Exponent -i/) as a function of (5/e at 
two-loop order for the RF, RB, and RP universality classes. 



model is organized in such a way that the single-replica 
p-loop diagrams correct the elasticity only to {p + 1) or- 
der. Since i?2(u) remains analytic along the FRG flow, 

we have R^\o) > 0, and therefore, -tp < 0. The corre- 
sponding values of the exponent ip computed for the RF, 
RB, and RP universality classes using the FPs found in 
Ref. m are shown in Fig. IH 



D. Roughness exponent to two- loop order 

We now show how one can calculate the roughness ex- 
ponent to second order in e and S knowing only the 
exponent ifj computed to second order in Sec. IIII CI To 
that end, we do not need the whole FRG to two-loop 
order. Let us start with the RB universality class. The 
roughness exponent is fixed by a stable RB FP solution 
of Eqs. (PO)) and (PT|) which decays exponentially fast for 
large u. The equations possess both the SR RB FP with 
R2{u) = and the LR RB FP with R2{u) ^ 0. The 
roughness exponent corresponding to the SR RB FP is 
known to second order in e and reads^^iSl 

CsRRB = 0.208298e + O.OOGSSSe^ + 0{e^). (31) 

Despite the smallness of the two-loop correction, the es- 
timation of the exponent in d = 1, Csrrb = 0.6866 given 
by Eq. (|3ip . visibly differs from the known exact result 
2/3. One can improve the accuracy of C by use the Pade 
approximant [2/1] involving also the unknown third order 
correction. Tuning the latter in order to reproduce the 
exact result 2/3 for e — 3, we end up with the expression 



*imp 



SRRB 



0.208298£ + 0.040017e2 
l + 0.159192e 



(32) 



which is expected to be fairly accurate for < e < 3. 

We now focus on the LR RB FP with i?2(u) ^ 0. 
We can integrate both sides of flow equation (pij) over 
u from to oo taking into account that for RB disorder 



a 2 



one-loop 
two-loop 
■ improved 



SRRBFP 



LRRBFP 



d 



FIG. 5: (Color online) Stability regions of SR and LR RB 
FPs on plane (d, a). The borderline between regions is given 
by equation ("lrrb = Csrrb (on the exact crossover line also 
•4^ — 0). The shown lines are computed using the one-loop, 
two-loop, and improved results. The circle is a point on the 
exact crossover line. 



R2{u) decays exponentially fast. Since for RB disorder 
the integral du i?2 (u) is nonzero, we can determine 
the roughness exponent Clrrb = S/5 to first order in 
e and 6. Fortunately, one can go beyond the one-loop 
approximation. Indeed, the direct inspection of diagrams 
contributing to the flow equation pT|) shows that the 
higher orders can only be linear in even derivatives of 
i?2 (m) . The only term which is linear in i?2 (u) comes from 
the renormalization of elasticity and can be rewritten as 
2ipR2{u) to all orders. Hence, we have to all orders 



((5-5C + 2V') / dui?2(u), (33) 
'o 



di / duR2{u) 



and as a consequence, duR2{u) is exactly preserved 
along the FRG flow resulting in the exact identity 



LRRB — 



5 + 2iP 



(34) 



Substituting Eq. ([50)1 into Eq. , we obtain the rough- 
ness exponent Clrrb to second order in e and S. Before 
we proceed to compute the exponents, let us to check 
stability of the SR and LR RB FPs. As was shown in 
Ref. [3^, the SR RB FP is unstable with respect to LR 
disorder if Clrrb > Csrrb. To one-loop order, this gives 
that the SR RB FP is stable for 6 < 1.0415£. Equating 
(p4)) and ((32|) . we can compute the stability regions to 
second order in e and S (see Fig. [5]). The alternative way 
to determine the crossover line relies on the requirement 
that the exponent ^ is a continuous function of e and 
S. It is zero in the region controlled by the SR FP, and 
therefore has to vanish when approaching the crossover 
line from the LR stability region. Since the ip is of sec- 
ond order in e and S, the ip criterion at two-loop order 
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FIG. 6: (Color online) RB disorder: roughness exponent as 
a function of 5 for e = 1 and e = 2. The solid lines are 
computed using the improved estimation (|35p and (|36p of the 
exponent ip. The dashed lines are computed using ip which is 
given by Eq. (|30p and shown in Fig. |4l The dotted line is the 
one-loop result (^lrrb = 5/5 which does not depend on e. 



gives the same stability regions as the roughness expo- 
nents equating at one-loop order. However, we can sig- 
nificantly improve the latter if we take into account that 
'0 = on the crossover line. The resulting crossover line 
is shown in Fig. [5] We can also improve the two-loop 
estimation of ■0. To that end, we write down a formal 
expansion of i/' in e. 



(35) 



The function fi{x) is basically the function shown in 
Fig. [4] We now tune the function f2{x) in order to make 
= on the crossover line and find 



0.159192a; - 0.200087 



X - 



1.04149 



(36) 



Using Eqs. ((35|) and ((36)) . we compute the roughness ex- 
ponent Clrrb as a function of 5 for e = 1 and e = 2 (see 
Fig. E]) . Unfortunately, the accuracy rapidly decays with 
e, so that estimation of the roughness exponent for e = 3 
is very difficult and postponed to Sec. |Vl 

Similar to the case of RB disorder, one can show that 
the roughness exponent at the LR RF FP is exactly given 

by 







LRRF 



5 + 2^; 



(37) 



and the crossover line between the SR and LR RF FPs 
is exactly given hy 5 = e. 

IV. RESPONSE TO TILT 

In this section, we study the response of a d- 
dimensional elastic object to a small tilting force tending 



to rotate the object in the plane {xi,u). The tilting force 
can be incorporated into the Hamiltonian as follows: 



nh[u] = lilu] - h / d'^xd.. 



(38) 



Such a force can be caused, for example, by a tilt of the 
applied field in superconductors or by tilted boundary 
conditions in the case of interfaces. For superconductors, 
we have h — (f)QH±/{4:'K), where H± is the component 
of the applied magnetic field transverse to the flux lines 
directed along xi and 0o is the magnetic flux quantum. 
Since we restrict our consideration to the case iV = 1, our 
results can be applied only to flux lines confined in (1-1-1) 
dimensions. However, the methods we use here can be 
extended to general N, and therefore applied to vortices 
in (2 -|- 1) dimensions. 

We focus on the response of the system to a small field 
h, which can be measured by the average angle between 
the perturbed and unperturbed orientations of the object 
in the {xi,u) plane: i!){h) := dx^Ux- In the absence of 
disorder the straightforward minimization of the Hamil- 
tonian leads to the linear response: ^{h) — h/c. To study 
the effect of disorder, it proves more convenient to work 
in the tilted frame: + -dxi. The corresponding 

Hamiltonian is 

f [l 

Hh[u] ^ j d'^x -'^c^idx^uxf + V[x,Ux + dxi) 



-{h - cid)dxiU3 



(39) 



where the field u satisfies dx^Ux = 0. Note that due 
to the violation of the STS symmetry, the tilted system 
can exhibit anisotropic effective elasticity even if the bare 
elasticity and disorder are isotropic. We now show by 
simple power counting that a finite tilt does introduce a 
new length scale in the problem which can be associated 
with the correlation length defined through the connected 
two point correlator. 



Q{x; ■&) = dxj^u{0)dxiu{x) ~ exp(-a;i/^^). 



(40) 



Indeed, upon scaling transformation x bx, u — > 
the arguments of the disorder term in Hamiltonian (j39p 
scale like V{bx, f'Ux+'&bxi). Comparing two terms of the 
last argument, we find that finite i? changes the character 
of disorder correlator above the length scale 



(41) 



diverging for ?9 provided that C < 1- Below one 
can neglect the tilt, while above the dependence on 
Ux is completely washed out and the i? term starts to 
suppress the correlation of disorder along xi. Thus, 
serves as the correlation length along xi, and therefore 
ci does not get renormalized beyond this scale. In the 
next two sections, we investigate the difference in the 
response to tilt for anisotropically distributed extended 
defects and isotropic LR-correlated disorder. 
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A. Response in the presence of columnar disorder 

Here, we extend the previous one-loop FRG 
studies^^i^ of elastic systems in the presence of columnar 
disorder to two-loop order and proceed to describe the 
transverse Meissner physics in a quantitative way. We 
consider the model with e^-dimensional extended defects 
introduced in Sec. Ill Al We take a = c\\ {i = 1, Sd) and 
we are also free to put Ci ^ 1 {i = Sd + l, ■■■,d) since they 
do not get corrected by disorder. Simple power counting 
shows that the upper critical dimension of the problem 
is due = 4 + e^. We use the dimensional regularization of 
integrals with a, e = A — d + £d expansion. The FRG flow 
equations to two-loop order read 

deRiu) = (e - 4:C)R{u) + (uR'iu) - TR'^^^u) 

+ ii?"(y)2 _ R"{0)R"(u) - ^R"'{0+)^R"{u) 

+ ^iR"{u)-R"{0))R"\uf, (42) 

a^lncii =R^^\0) + R^^\0f + 2R'"{0)R'-^\0), (43) 

9flnr=-0T-yi?<^^(O) + O(ii''), (44) 

deh = c;/'Aoe-^ii""(0+) + 0(i?2). (45) 

where 9t — d — 2 + and Aq is the bare cutoff. In 
Eq. ([^5]) , h is the coefficient in front of the term 

[ d^x\V\\ul\ (46) 

a '' 

ultimately generated in the effective Hamiltonian along 
the FRG flowi^ The correction to h is strongly UV di- 
verging, and thus is nonuniversal. Note that the flow 
equation for the generalized columnar disorder (|42p coin- 
cides to all orders with that for pointlike disorder up to 
change e ^ e. 

Let us start from the analysis at T = 0. The flow 
picture very resembles that for the depinning transition 
with tilt "i?, longitudinal elasticity C|| , and tilting field h 
playing the roles of velocity, friction, and driving force, 
respectively. For d < 4, + ed^ the running disorder corre- 
lator R^^^ (u) blows up at the Larkin scale 

Ic^ l\n{l + ). (47) 

Consequently, the longitudinal elasticity diverges at zero 
tilt z9 in a way similar to mobility divergence at the 
depinning transition in quasistatic limit. Beyond the 
Larkin scale I > Ic, the i?(^)(u) develops a cusp at origin, 
i?"'(0+) > 0. The term (gS]) is generated in the effective 
Hamiltonian and the i?^'*-'(0) changes its sign from posi- 
tive to negative. The latter leads to a power law decay 
of the longitudinal elasticity cy cx L^'^' with 

^ = -i?*(4)(0) - i?*W(o)2 _ 2i?*"'(0)i?*(5)(0), (48) 



z9 




FIG. 7: (Color online) Schematic plot of the elastic object 
response to a transverse field for uncorrelated disorder (dot- 
ted line), columnar disorder (dashed line), and LR-correlated 
disorder (solid line). 

where R*{u) is a FP solution of the flow equation (|42p . 
Similar to the threshold force generation at the depin- 
ning transition, term (j46p reduces the tilting force and 
generates the critical tilting force he- The flow equation 
PS)) allows us to estimate the nonuniversal value of he- 
Integrating Eq. (|45| up to large scales, we find 

he = ° ^l^^^ R*"'{0+). (49) 

We now in a position to compute the exponent which 
we define as 

i?cx (/i- /ic)"*. (50) 

To that end, we renormalize the equilibrium balance 
equation h ~ he — ci{L)'d up to the scale L = ^,9 at 
which the elasticity ci stops to get renormalized. Using 
Eq. (|iT1) . we obtain the exact scaling relation 

The exponents ip, and computed to second order 
in s for different universality classes are summarized in 
Tab. D Note that expansions in s are expected to be 
Borel nonsummable, and thus ill behaved for high orders 
and large e. In this light, the using of exact relation 
(jSip may be more favorable than the expansions given 
in the last column of Table HI Systems described by the 
RP universality class exhibit slow logarithmic growth of 
displacements 

{ux^ - ua)'^ = Adlnxj_, (52) 

where u^^ := L~^'' J d'^''x\\ Ux- The universal amplitude 
can be easily deduced from the results for the uncorre- 
lated disorder and to two loop reads 

= ^ + + (53) 
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TABLE I: Critical exponents for elastic systems with gener- 
alized columnar disorder computed to two-loop order. 





c 




<^ 




RP 





3^ ^ 9^ 


1 




RF 




9^ ' 162'^ 


1 


9^ 18^ 


RB 


0.208298e 


0. 2639025 


1 


- 0. 2639025 




+ 0.006858e^ 


+ 0.053615£^ 




0.038941 



Here, we restore the dependence on the flux line den- 
sity riQ fixing the period of R{u) to I/uq. In contrast to 
the Bragg-glass, the weak Bose-glass survives in d = 2. 
Indeed, for uncorrelated disorder in d = 2, the tem- 
perature turns out to be marginally relevant, so that 
the system has a line of FPs describing a super-rough 
phas e with ano malous growth of the two-point correla- 
tion {ux — uq)^ = A{T) In^ X + O(lnx).'*^ According to 
Eqs. (j55p and (|56p for columnar disorder the tempera- 
ture vanishes at finite, though a very large scale 



where we have fixed the period to 1. The logarithmic 
growth of displacements corresponds to a slow power- 
law decay of the translation order, and thus should lead 
to Bragg peaks unexpected for a strongly pinned Bose- 
glass. The system under consideration shares features of 
the Bragg glass, such as power-law decay of the trans- 
lation order, and the strong Bose glass, namely, the di- 
verging tilt modulus (transverse Meissner effect). One 
can expect this behavior for a weak Bose-glass which is 
pinned collectively. Recently, such a glassy phase called 
Bragg-Bose glass was observed in numerical simulations 
of vortices in bulk superconductor at low concentrations 
of columnar disorder and low temperatures 4^:^ 

A finite temperature T > rounds the cusp of the 
running disorder correlator Ri{u), so that in the bound- 
ary layer u ~ T;, it significantly deviates form the FP 
solution and obeys the following scaling formi^ 



i?r(«) = i?r(o)-Tz 1 - vi + {ux/Ti) 



(54) 



where x — \R*"'{0)\- However, as was pointed in Ref. 33 
the flow equations for columnar disorder have a remark- 
able feature in comparison with uncorrelated disorder. 
Indeed, substituting the boundary layer scaling (j54p in 
the temperature flow equation (|44p . we obtain 



(55) 



As follows from Eq. ([55]), the effective temperature T; 
vanishes at a finite length scale Lioc = e'i°'^/Ao, 



loc 



9t 



In 1 



29tTq 



(56) 



so that the localization effects are settled only on scales 
larger than ^loc > Ic- 



B. Interacting disordered bosons in (1 + 1) 
dimensions 

Let us discuss the special case of flux lines in (1+1) di- 
mensions which the qualitative phase diagram is shown 
in Fig. [T] The transverse Meissner physics for collectively 
pinned weak Bose glass and small tilt angles d = Bj_/ B\\ 
can be explored using the results obtained in the previous 
section for the RP universality class with e = 4—2+1 = 3. 



Hoc 



2To 



(57) 



Unfortunately, the large value of e makes estimation of (p 
extremely unreliable. Indeed, the expansion in e shown 
in Table U leads to a zero value of 0. The exact scaling 
relation (jSip with -0 computed using the expression from 
Table |T] gives 



(1 + 1) 



- (one loop), 



(1 + 1) = - (two loop). 
3 



(58) 



The one-loop result reproduces the estimation (j) — 1/2 
given by heuristic random walk arguments based on the 
entropy of flux lines wandering in the presence of ther- 
mal fluctuations.^^ The model of vortices wandering in a 
random array of columnar defects can be mapped onto a 
quantum problem of disordered bosons.^ One can regard 
each vortex as an imaginary time world line of a boson, 
so that the columnar pins parallel to vortices become 
quenched pointlike disorder in the quantum problem. 
The transverse magnetic field H± will play the role of 
an imaginary vector potential h for the bosons^^ so that 
the bosonic Hamiltonian turns out to be non-Hermitian: 



n 



2m 
1 



dxip^ (x) {dx ~ h)'^ il;{x) + / dxV{x)h{x) 



+ — j dxdx' n{x)U{x ~ x')h{x'). 



(59) 



Here, ili^x), ipi^) tli6 bosonic creation and annihi- 
lation operators, and n{x) = ^^{x)'ijj{x) is the density 
operator. U (x) is a short-range repulsive interaction po- 
tential between bosons with the strength Uq = J dx U{x). 
The disorder is described by a time-independ ent G aus- 
sian random potential V{x) with zero mean V{x) — 
and short-range correlations V{x)V{x') = Vo6{x — x'). 
We can pass to a quantum hydrodynamic formulation 
of model ([55)1 expressing everything though the bosonic 
fields 0{x) and ^p{x) which satisfy the canonical commu- 
tation relation^ 



[dMx),0{x')]^~i5{x-x'). 



(60) 
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For bosons with average density no, this gives 



TT g 



/oo 
dxv{x) 



(61) 



p— — oo 



where in the disorder part we have retained only the lead- 
ing contributions coining from the backscattering on im- 
purities. The forward scattering term can be eliminated 
by a shift of the phonon field 9{x) which does not depend 
on the time t, and thus, this term does not contribute to 
the current J ~ dt9{x). The Luttinger liquid parameter 
g and the phonon velocity Vp are given by 



mllQ ' 



UoUo 



(62) 



The imaginary time (r = it) action corresponding to 
Hamiltonian (j6ip for a particular distribution of disorder 
can be derived using the canonical transformation 



dxdT[n ~ ine{x)drO{x)]. 



(63) 



Here, drOix) = {i/h)[H,9{x)] = -igcng{x)/{TTh) and 
Hg = hdxf{x) is the momentum conjugate to 9{x) which 
is given by Eg. ([50]) . Averaging 1^ Qygj. disorder by 
means of the replica trick and keeping only the most rel- 
evant terms, we obtain the replicated action 



S 



y dxdr 



2^ 



-{drea{x,T)f 



+ Vp{dx9a{x,T)Y 



hd- 



r9a{x,T)] 



^ / dxdrdr' cos{2TT[9a{x,T) ~ 9b{x,T')]} . (64) 



The imaginary time action (|64p is identical to the Hamil- 
tonian of periodic elastic system with columnar disorder 
^ . The imaginary time plays the role of the longitudi- 
nal coordinate r < — > x\\ which is parallel to columnar 
pins. The Planck constant stands for the temperature 
h < — > T, and the phonons are related to the dimension- 
less displacements field 9{x,t) = —nou{x). There is the 
following correspondence between quantities in the vor- 
tices and bosons problems^ 



(65) 



The vortex tilt angle d caused by the transverse field 
H±_ corresponds to the boson current J = {—i)d'H/dh 
induced by the imaginary vector potential h. For h = 0, 
the disordered bosons undergoes a superfluid-insulator 
transition at g = 3/2. This determines the temperature 
Tbg, such that ^(Tbg) — 3/2, above which vortices form 
a liquid (see Fig. [1]) . It is known that in one dimension 



there is no difference between bosons and fermions, and 
both types of particles are described by the Luttinger 
liquid (|6ip . In particular, the hard-core bosons can be 
mapped onto free fermions that corresponds to a special 
value of the Luttinger parameter g{T*) = 1, which de- 
fines the temperature T*. In Ref. 14, the mapping onto 
free fermions was used to study the transverse Meissner 
effect in (1-t-l) dimensions. The free fermions on a lattice 
is described by the tight-binding model, 



n = 



|^-Wi(c|ci+ie -I- cl+jCic'') -f {ti - /i)cjci 



(66) 



where c^^, c are on site fermion creation and annihilation 
operators, and /i is the chemical potential. Wi is a ran- 
dom hopping matrix element and is a random pinning 
energy. In Ref. both cases, the random pinning and 
the random hopping models, were studied using the exact 
results for the Lloyd model and the strong-randomness 
real-space RG, respectively. It was found in both cases 
that J ^ h — /ic, i.e., </) = 1, that significantly differs 
from the FRG prediction ([55)1 . The difference can be at- 
tributed to that the free fermions analog is limited to a 
special point g = 1 (T = T*), while the FRG prediction 
may be valid only for low temperatures since it is con- 
trolled by the zero-temperature fixed point. The corre- 
spondence between the temperature and the Planck con- 
stant in both problems reflects that the zero-temperature 
FRG FP may have a counterpart in the quantum problem 
in the form of an instanton solution. This may account 
for the consistency of the exponent computed by FRG 
and estimated using heuristic arguments of kink statis- 
tics. 

The high- Tc superconductor films grown by deposition 
often exhibit larger critical currents than their bulk coun- 
terparts due to the formation of dislocations running par- 
allel to the crystalline axis, and thus, they are natural 
candidates to verify the above results. However, as was 
discussed in Ref. [50, the picture may be more involved 
since the dislocation lines can meander or they can be of 
relatively short length that breaks up the Bose glass into 
pieces along the direction of the crystalline axis. 



Response in the presence of long-range 
correlated disorder 



We now consider the response to tilt in the presence 
of isotropic LR-correlated disorder. In contrast to the 
case of generalized columnar disorder the h term is not 
generated due to the analyticity of the LR part R2{u) 
of disorder correlator. Moreover, the elasticity remains 
finite along the FRG flow though it grows as a power 
law c ^ with ip < given by Eq. pO]) and shown in 
Fig. 21 As a consequence, there is no threshold transverse 
field: the systems is tilted for any finite tilting force. 
Renormalizing the balance equation ft = cii? up to the 
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scale given by Eq. pTj) . we see that the response to 
the tilting force h is given by a power law 



(67) 



with the exponent (j) > 1 defined by Eq. ([5T|l . The re- 
sponse to tilt in systems with uncorrelated, columnar 
and LR-correlated disorders is shown in Fig. [T] As one 
can see from the figure, the response of systems with 
LR-correlated disorder interpolates between the response 
of systems with uncorrelated and columnar disorder. In 
particular, we argue that in the presence of LR-correlated 
disorder, vortices can form a new vortex glass phase 
which exhibits Bragg peaks and vanishing linear tilt mod- 
ulus without transverse Meissner effect. We will refer to 
this phase as the strong Bragg glass. 

In analogy with the Bose glass, one can attempt to 
map the system with linear defects of random orientation 
corresponding to LR-correlated disorder with a — d ~ 1 
to a quantum system consisting of interacting bosons and 
heavy particles moving with random quenched velocities 
according to classical mechanics. 



V. KARDAR-PARISI-ZHANG EQUATION 
WITH TEMPORALLY CORRELATED NOISE 

In this section, we address the relevance of our results 
to the Kardar-Parisi-Zhang (KPZ) equation (and closely 
related Burgers equation), which describes the dynam- 
ics of a stochastically growing interface.^i The latter is 
characterized by a height function h{x, t), x which 
obeys the nonlinear stochastic equation of motion 



dth = vV'^h + ^(V/i)2 + -nix, t). 



(68) 



The first term in Eq. (|68p represents the surface tension, 
while the second term describes tendency of the surface 
locally grow to normal itself. The stochastic noise rj{x,t) 
is usually assumed to be Gaussian with short-range corre- 
lations. Here, we consider the noise with long-range cor- 
relations in both time and space. It is defined in Fourier 

{r]{k,uj)T]{k',uj')) = 2D{k,uj)5'^' {k + k')5{u + J), (69) 

with the noise spectral density function having power-law 
singularities of the form 



D{k,Lu)=Do + Dgk 



-2P„ 



-29 



(70) 



Such temporal correlations can originate from impurities 
which do not diffuse and impede the growth of the in- 
terface, while the space correlations can be due to the 
presence of extended defects. Since there is no intrinsic 
length scale in the problem, asymptotics of various cor- 
relation functions are given by simple power laws. For 
instance, the height-height correlation function scales like 



{[hix,t)-hix',t')f)^\x-x'\'^f 



\t~t'\ 

\x — x'p 



(71) 



where x is the roughness exponent and z is the dynamic 
exponent which describes the scaling of the relaxation 
time with length (do not mix it with the dynamic expo- 
nent z at the depinning transition, which is not used in 
this paper). 

Medina et al.^^ studied the KPZ equation with the 
noise spectrum (I70|) using the dynamical renormalization 
group (DRG) approach and here we adopt the notation 
introduced in their work. Let us briefly outline the results 
obtained in Ref. [s^ restricting ourselves mainly to the 
case d! — 1. The flow equations expressed in terms of 
dimensionless parameters = KdiX^D^/v^ and Ug — 
Kd'X^Dg/i'^ to one loop order read 

delniy = z - 2 + ^ + ^{1 + 2p){l + 26) sec(7r6'), (72) 

di\iiX = x + z-2 + UeOil + 2d) scc{t:9), (73) 
diUe^Ug[z{l + 20)-2x-l + 2p], (74) 

5,t/o ^ U„{z - 2x - I) + ^ + ^{1 + 40) sec(2^0) 
1, 



+ 2'^oC/9(l + 26l)sec(7r6 



(75) 



Note that the DRG calculations are uncontrolled, in the 
sense that there is no small parameter. For white noise 
(9 = 0), the KPZ equation is invariant under tilting of 
the surface by a small angle. The STS symmetry implies 
that the vertex A does not get corrected by the noise to 
all orders. This results in the exact identity 



XSR + ZSR 



(76) 



Besides the known SR FP with Ug = 0, the flow equa- 
tions ([7^ - (|75p are expected to have a different LR FP 
with Ug ^ Q. It was argued that the term Ug in the noise 
spectrum -D(fc, uj) acquires no fluctuation corrections: the 
scaling of Ug is completely determined by its bare dimen- 
sion so that Eq. (|7^ is exact to all ordersi^ This allows 
one to compute the exact critical value 9c = 1/6 (for 
/9 = 0) at which there is a crossover from the SR FP to 
the LR FP. For arbitrary p, the crossover to the LR FP 
happens at 



W + Ap> 1. 



(77) 



The term Ug becomes relevant and as follows from 
Eq. ((74|) the exact relation 



zlr(1 + 29) - 2xLR + 2p^l 



(78) 



holds at the LR FP. Let us for the moment ignore the 
noise correction to A in Eq. (j73p . This approximation 
restoring the STS is valid only for small 9 and yields 



z*{e,p) 



^ l + 49 + 2p 
3 + 261 ' 
1 + 46* + 2p 
3 + 261 ■ 



(79) 
(80) 



For large 9, one can expect a significant deviation of ex- 
ponents z and X from z* and x* ■ To gain insight into the 
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problem the authors of Ref. 'sS solved the flow equations 
([7^ - (|75p for finite and p = numerically. They found 
that the physical LR FP exists only for 6* < 1/4, while 
nothing special is physically expected at 9 — 1/4. It 
was argued that the problem is originated from infrared 
divergences of integrals and that infinite number of addi- 
tional terms generated in the noise spectral density under 
DRG: 



(81) 



Keeping track of renormalization of all Z?„, the authors 
of Ref. [S^ solved the truncated system of flow equations 
numerically and found that the critical exponents for p — 
can be fitted to 



Z 1.50-lr 



one-loop FRG/DRG 
■ two-loop FRG 
modified DRG 
self-consistent 
simulations 



0.2 0.3 





x{0) = 1.696* + 0.22, 
2x{e) + 1 



z{e) 



1 + 20 



(82) 
(83) 



We now revise the problem in the light of what has 
been learned in the previous sections. Using the well- 
known Cole-Hopf transformation Z — exp[(X/2i>)h] one 
can eliminate the nonlinear term in Eq. (|68p and obtain 
a diffusion equation in time-dependent random potential 

dtZ(x, t) = vV'^Zlx, t) + —-qlx, t)Z{x, t). (84) 

2iy 



The solution of Eq. (|84|) can be regarded as the partition 
function of a directed polymer (DP) of length t in {d' + 1) 
dimensions with ends fixed at (0,0) and {x,t): 



Z{x,t) 



x{t)=x r ^ 

Dx(t) exp<^ 
x(o)=o L ^ 



dt 



+7j{x{t),t)\} 



(85) 



with u = T/2c and A = 1/c. The DP is a one-dimensional 
{d — l,e = 3) elastic object with d' = N - dimensional 
target space. Thus, the time-dependent noise 'q{x, t) in 
the KPZ equation is mapped to the quenched disorder 
V in the DP picture. This gives the exact relation be- 
tween the dynamic exponent of KPZ problem and the 
DP roughness exponent which reads 



z(d') = l/C(d=l,iV = d')- 



(86) 



Spatial correlations in •q{x,t) corresponds to correlations 
of quenched disorder V in the directions transverse to the 
DP. As the exponent p varies from to 1, the quenched 
disorder interpolates between RB and RF universality 
classes. For example, the exponent z changes from 3/2 
to 1 for d' — 1 and white random noise [6 — 0). The sta- 
bility criterion assures that the LR FP in the FRG picture 
is stable if Clr > Csr- This implies that the noise tem- 
poral correlations in surface growth problem are relevant 
only if the corresponding dynamic exponents fulfill the 
condition zlr < ^SR- Note that this criterion is purely 



FIG. 8: (Color online) Dynamic exponent z for the KPZ equa- 
tion [d! — 1) with temporally correlated noise (p — 0) com- 
puted using different technique: one- loop FRG/DRG is given 
by Eq. ([79}; two-loop FRG is given by Eq. ^} (note that 
large value £ = 3 does not allow for accurate computation 
of tp and z at large 6); modified DRG is given by Eq. (I83|l : 
self-consistent approximation of Ref. [sj and simulations of 
Ref [13. 



based on the mapping between the DP and KPZ prob- 
lems. Since zsR{d' = 1, p = Q) = 3/2 the exponent ([55]) 
computed using the modified DRG violates the criterion 
of the LR FP stability, and thus is ruled out. Substi- 
tuting the roughness exponents computed using FRG for 
the RB {p — 0) and RF {p — 1) universality classes into 
Eq. ([55]) and relating (5 = 3-1- 20, we obtain the exact (for 
p = 0, 1 and presumably for any p) identity 



^LR 



2p 



3 + 20 + 11; 



(87) 



To one-loop order in FRG, i.e., for ip = 0, exponent 
([87| coincides with the estimation given by DRG ([79]) for 
small 0. Though the exponent ip has been computed in 
Sec.|niC]for the RB {p = 0) and RF (p = 1) universality 
classes to two-loop order in a controllable way, the large 
value £ = 3 of the expansion parameter describing the 
DP problem makes the estimation of ip highly unreliable. 
Nevertheless, since ^ = is zero on the crossover line 
between the LR and SR FPs, we can determine this line 
exactly for p = 0,1 from equation zlr < zgR =3/2 that 
leads back to Eq. ([77]) . Taking into account that is 
nonpositive for columnar disorder, we obtain the lower 
and upper bounds on z (0) for p = and 6* e [ii] as 



(88) 



The critical exponent z computed using FRG, DRG, and 
measured in numerical simulations of Ref. is shown in 
Fig. [H The KPZ equation with temporally correlated 
noise was also studied using a self-consistent approxi- 
mation (SCA).^'* The SCA equations have two strong- 
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coupling solutions. The first one exhibits a crossover- 
like behavior a,t 9 = ^ and corresponds to the one-loop 
FRG prediction. The second solution, which is consid- 
ered to be dominant, leads to a smooth dependence of 
z on. 9 shown in Fig. [S] Both the SCA solutions are in 
agreement with the FRG prediction that the exponent 
2 is a decreasing function of 9, while the modified DRG 
suggests that z increases with 9. However, the second 
SCA solution considered to be dominant does not satisfy 
bounds and thus is ruled out. 

Let us generalize identity ([75]) to the case of tempo- 
rally correlated noise. Note that the solution of the KPZ 
equation h{x,t) gives the free energy of DP ([85]) . The 
free energy per unit length /(i?) of the DP tilted by the 
transverse field to the angle d can be written as 
= /(O) + ct9" — i?_L The naive elastic approxi- 
mation suggests a = 2. In order to take into account the 
renormalization of elasticity, we determine the exponent 
a from the condition that at equilibrium the response to 
the field H± is d H'^. This fixes a = l-f 1/0 with 
given by Eq. ([51]) . Then the total free energy of the DP 
of length t can be written as a function of the free end 
coordinate x as follows: 

h{x) ^tf{0) + tc{x/t)" -x'^. (89) 

The last term in Eq. ([89|l describes the typical fluctuation 
of the free energy due to the disorder and is given by 
Eq. ([71]) . Balancing the last two terms of Eq. ([89]) and 
using Eq. ([51]) . we obtain the exact scaling relation 

X + z^2-i;z, (90) 

which holds at the LR FP as well as at the SR FP. At 
the SR FP V = 0, so that Eq. ^ reduces to the STS 
identity Excluding ip from Eqs. ^ and ((87|), we 

arrive at the relation ([75]) valid at the LR FP. 



glass which is expected in type-II superconductors with 
columnar disorder at small temperatures and at high vor- 
tex density which exceeds the density of columnar pins. 
The weak Bose glass is pinned collectively and shares 
features of the Bragg glass, such as a power-law decay of 
translational order, and features of the strong Bose glass, 
such as a transverse Meissner effect. For isotropic LR- 
correlated disorder, the linear tilt modulus vanishes at 
small fields leading to a power-law response I? ~ /i"^ with 
(j) > 1. The response of systems with LR-correlated dis- 
order interpolates between the response of systems with 
uncorrelated and columnar disorder. We argued that in 
the presence of LR-correlated disorder vortices can form 
a strong Bragg glass which exhibits Bragg peaks and a 
vanishing linear tilt modulus without transverse Meissner 
effect. The elastic one-dimensional interface, i.e., the di- 
rected polymer, in the presence of LR-correlated disorder 
can be mapped to the Kardar-Parisi-Zhang equation with 
temporally correlated noise. Using this mapping, we have 
computed the critical exponents describing the surface 
growth and compared with the exponents obtained using 
dynamical renormalization group, self-consistent approx- 
imation, and numerical simulations. 
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VI. SUMMARY 

We have studied the large-scale behavior of elastic sys- 
tems such as interfaces and lattices pinned by correlated 
disorder using the functional renormalization group. We 
consider two types of disorder correlations: columnar dis- 
order generalized to extended defects and LR-correlated 
disorder. Both types of disorder correlations can be pro- 
duced in real systems, for example, by subjecting them 
to either static or rotating ion beam irradiation. We 
have computed the critical exponents to second order in 
e = 4 — d and 5 = 4 — a for LR-correlated disorder and 
to second order in e = 4 — e -I- for ed-dimensional ex- 
tended defects. The correlation of disorder violates the 
statistical tilt symmetry and results in a highly nonlinear 
response to a tilt. In the presence of generalized colum- 
nar disorder, elastic systems exhibit a transverse Meiss- 
ner effect: disorder generates the critical field he below 
which there is no response to a tilt and above which the 
tilt angle behaves as i) ^ {h — h^)'^ with a universal ex- 
ponent (j) < 1. The periodic case describes a weak Bose 



APPENDIX A: CORRECTION TO ELASTICITY: 
TWO-LOOP DIAGRAMS 

In this Appendix, we calculate diagrams shown in 
Fig.|3]keeping only the terms which correct the elasticity. 
Here, we set gi{x) := 5'^{x) and g2{x) := g{x). Diagram 
a yields 

[o\aP = ~9T / RqcWx ~ ^x-yi-y^-ys) 

•^yi, J/2, J/3 

xRf\u^^y^ - Ux^y^^y^)gaiyi +y2+ ya) 

3 

^9piy2)Y[C{yih (Al) 

1=1 

where a, P = 1,2 and C{x) is given by Eq. ©. Using 
the short distance expansion (|23p . we obtain 

3 

HS = -74;(V«.)'i?i'H0)<(0) / flCim) 
^"■^ ''{y}i=i 

X {yi + y2 + y^fgaiyi + y2 + yz)gfi{y2) (A2) 
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and 



Straightforward analysis shows that 



3 



l(2) ^ 

J a/3 4^j.V 



(A3) 



Note that the term R'^'{0)R'^'{0) does not contribute 
since LR disorder i?2(u) remains an analytic function 
along the FRG flow, while all diagrams with a = /3 — 1 
are zero due to the STS. We will neglect similar terms in 
what follows. For diagram b, we have 



[b] 



af3 



2T 



xRf\ux ~ Ux_yjga{y2 - yi) 



{y3)U^iyi 



(A4) 



1=1 



Applying the short distance expansion (|23p . we arrive at 

3 

X (yi + y2f9a{yi + y2)gp{y3) (A5) 

and 

3 

■^{v}i=i 

xyl9a{yi + y2 )gp{.y?,)- ( A6) 

Diagram [c] gives 



[Cja/3 = 



^x—yi '^x~yi—y2—y3 )9a{yi +y2) 
3 

^9(3{y2 + y3)Y[C{yi). (A7) 

1=1 

After short distance expansion, we find that diagrams [c] 
give rise to elasticity correction only for a = f3 = 2 which 
reads 

3 

[cU = -^{ynxrRi'H0)Rf\0) I t{C{yi) 
X [{yi + y2) ■ (y2 + y3)]9a{yi + y2)9i3{y2 + ys)- (A8) 



i(i) 

J la 
l(l) 



l(2) 



[^]i'j = [d^=0, (a = 1,2), 



l(2) 



L"'J21 [^]21^ — 0, 
[C]ll = [C]i2 = [c]21 =0. (A9) 

We now compute the integrals combining them in pairs: 



I 



(gi + 92)2 + I 



a—d—2 a—d 

.2 , ?.2 . (AlO) 



(9? + !)^ 



where we have included l/c^ in redefinition of Ri (u) . The 
integral over qi in Eq. (|A10p is of order 0(l/e) so that 



[a] 22 + [^]22 is finite and does not correct elasticity at 
two-loop order. Other diagrams give 



l(2) 



+ [b] 



(2) _ c(Vu^)2 (g - d){a - 2) ^^^^-^ 



2T 2d 

x(Ji- J2)i?'i'(0)i?f^(0), 



(All) 



l(2) 



[6] (2) ^ c(Vu^)2 (a_rf)(a_2) 



2T 



2d 



:(J3- J4)i?2(0)i?f (0) 



(A12) 



where we have defined the following two-loop integrals: 



Ji 

J2 



a-d-2 
I2 



Kl 



{l + qlY[l + {qi+q2Y] 5{5 + s) 

a-d-2 
H2 



{l + qly{l + ql) 
{l + qj)2[ll{q,+q,)2] 26^ 



= §2+0{e-\5-^) 



a—d a—d— 2 



,(1+^2)2(1+^2) 

The last diagram 

(a-d)2c (4) 2 



(A13) 



H22 



2dTm 

[<?! • 92] qt 



a-d-2 a-d-2 
I2 



M [iqi+q2r + lKql + lM + l) 



(A14) 



is finite in the limit e,S —^ 0, and thus does not correct 
elasticity at two-loop order. 
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